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Abstract 

The experimental data on net proton and net A-hyperon spectra obtained by the 
NA35 Collaboration, as well as the inclusive densities of A and A obtained by NA49, 
NA57, and STAR collaborations, are compared with the predictions of the Quark- 
Gluon String Model. The contributions of String Junction diffusion, interactions with 
nuclear clusters, and the inelastic screening corrections are accounted for. The level 
of numerical agreement of the calculations with the experimental data is of about 
20—30%. The predictions for LHC are also presented. 



PACS. 25.75.Dw Particle and resonance production 
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1 Introduction 



The Quark-Gluon String Model (QGSM) [TJ is based on the Dual Topological linea- 
rization, Regge phenomenology, and nonperturbative notions of QCD. This model is 
successfully used for the description of multiparticle production processes in hadron- 
hadron [2j [3j HJ [5] , hadron- nucleus [6j [7], and nucleus- nucleus [H HI [10] collisions. In 
particular, the rapidity dependence of inclusive densities of different secondaries (n ± , 
K , p, and p) produced in Pb-Pb collisions at 158 GeV/c per nucleon were reasonably 
described in ref. [10]. In the present paper we consider the yields of p and p, as well 
as A and A, produced in the collisions of different nuclei at CERN SpS and RHIC 
energies. 

In the QGSM high energy interactions are considered as proceeding via the exchange 
of one or several Pomerons, and all elastic and inelastic processes result from cutting 
through or between Pomerons [TTj . Inclusive spectra of hadrons are related to the 
corresponding fragmentation functions of quarks and diquarks, which are constructed 
using the Reggeon counting rules [T2] . 

In the case of interaction with nuclear target, the Multiple Scattering Theory 
(Gribov-Glauber Theory) is used. It allows to consider the interaction with nucleus as 
the superposition of interactions with different numbers of target nucleons [131 EH HS1 

us]. 

Also in the case of nucleus-nucleus collisions the Multiple Scattering Theory allows 
to consider the interaction as the superposition of separate nucleon-nucleon interac- 
tions. Though in this case the analytical summation of all the diagrams is impos- 
sible |17j . the signiificant classes of diagrams can be analytically summed up in the 
so-called rigid target approximation [IS] which is used in the present paper. 

The significant differences in the yields of baryons and antibaryons in the central 
(midrapidity) region are present even at high energies. This effect can be explained 
[H El HSl [201 [2TJ [221 [231 121] in QGSM by the special structure of baryons consisting of 
three valence quarks together with a special configuration of gluon field, called String 
Junction [251 1251 (23 121] ■ 

One additional contribution comes from the coherent interaction of a projectile with 
multiquark clusters inside the nuclei. The existance of these interactions is confirmed 
by the presence, with not such a small probability, of a cumulative effect [25] . 

These contributions were incorporated into the QGSM in [301 [31], an d they allow 
to describe a number of experimental facts. 

At very high energies the contribution of the enhancement Reggeon diagrams be- 
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comes important, leading to a new phenomenological effect, the supression of the in- 
clusive density of secondaries [32] in the central (midrapidity) region. 

In this paper we present the description of p, p, A, and A production on nuclear 
targets at CERN SpS and RHIC energies. 

2 Baryon/antibaryon asymmetry in the QGSM 
2.1 General approach 

The QGSM [H [2J [3] allows one to make quantitative predictions for different features 
of multiparticle production, in particular, for the inclusive densities of different secon- 
daries, both in the central and in the beam fragmentation regions. 

In QGSM, each exchanged Pomeron corresponds to a cylindrical diagram, and thus, 
when cutting one Pomeron, two showers of secondaries are produced (see Fig. 1 a,b). 

The inclusive spectrum of a secondary hadron h is then determined by the con- 
volution of the diquark, valence quark, and sea quark distributions, u(x,n), in the 
incident particles, with the fragmentation functions, G (z), of quarks and diquarks 
into the secondary hadron h. Both the distributions and the fragmentation functions 
are constructed using the Reggeon counting rules. 

In particular, in the case of n > 1, i.e. in the case of multipomeron exchange, the 
distributions of valence quarks and diquarks are softened due to the appearance of a 
sea quark contribution. There is some freedom [B] in how to account for this effect. 
In principle, the valence and sea quarks can depend on n in a different manner, for 
example: 
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where every distribution Ui(x,n) is normalized to unity, d is a parameter, and the 
values of m 1 and m 2 can be found from the conditions: 

(x uu ) + (xd) + 2(n - 1) ■ (x s ) = 1 
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Figure 1: (a) Cylindrical diagram representing a Pomeron exchange within the DTU classification 
(quarks are shown by solid lines); (b) One cut of the cylindrical diagram corresponding to the single- 
Pomeron exchange contribution in inelastic pp scattering; (c) One of the diagrams for the inelastic 
interaction of one incident proton with two target nucleons N\ and N2 in a pA collision. 



(xud) + <ar«) + 2(n - 1) • (x s ) = 1. (2) 

The details of the model are presented in [TJ [2j [3j [19]. The averaged number of 
exchanged Pomerons (n) pp slowly increase with the energy. The Pomeron parameters 
have been taken from [3]. 

For a nucleon target, the inclusive rapidity, y, or Feynman-x, Xp, spectrum of a 
secondary hadron h has the form [TJ: 

^ = — 'T- =tw n .^ n (x)+w D .^), (3) 

where the functions determine the contribution of diagrams with n cut Pomerons, 

w n is the relative weight of this diagram, and the term wo ■ (j)%(x) accounts for the 
contribution of diffraction dissociation processes. 



In the case of pp collisions: 



= /*(x + ,n)./J(x_,n)+/J(x + ,n)./i(x_,n)+2(n-l)./*(x + ,n)./*(a;-,n) , (4) 

z± = ^[^/4rr4/s + a: 2 ±a;] , (5) 

where f qq , f q , and f s correspond to the contributions of diquarks, valence quarks, and 
sea quarks, respectively. 

These contributions are determined by the convolution of the diquark and quark 
distributions with the fragmentation functions, e.g., 

fq(x+,n) = / u q (x u n) ■ G h q (x + /xi)dxi . (6) 

In the calculation of the inclusive spectra of secondaries produced in pA collisions 
we should consider the possibility of one or several Pomeron cuts in each of the v blobs 
of proton-nucleon inelastic interactions. For example, in Fig. lc it is shown one of the 
diagrams contributing to the inelastic interaction of a beam proton with two target 
nucleons. In the blob of the proton-nucleonl interaction one Pomeron is cut, while in 
the blob of the proton- nucleon2 interaction two Pomerons are cut. The contribution 
of the diagram in Fig. lc to the inclusive spectrum is 

= 2.W M (2).^.<».{/*(x + ,3)./J(x_,l) + 

® 'prod F 

+ f h q {x + , 3) • 1) + 3) • [/* 2) + 2) + 

+ 2 ■/*(*_, 2)]} , (7) 

where W p a{2) is the probability of interaction with namely two target nucleons. 

It is essential to take into account all digrams with every possible Pomeron config- 
uration and its corresponding permutations. The diquark and quark distributions and 
the fragmentation functions are the same as in the case of pN interaction. 

The total number of exchanged Pomerons becomes as large as 

(n) p A ~ (v) p a ■ (n) P N , (8) 

where (v) p a is the average number of inelastic collisions inside the nucleus (about 4 for 
heavy nuclei at SpS energies). 



5 



The process shown in Fig. lc satisfies [131 [HJ [T5J UB] the condition that the absorp- 
tive parts of the hadron-nucleus amplitude are determined by the combination of the 
absorptive parts of the hadron-nucleon amplitudes. 

In the case of a nucleus-nucleus collision, in the fragmentation region of projectile 
we use the approach [HO [10], where the beam of independent nucleons of the projectile 
interact with the target nucleus, what corresponds to the rigid target approximation 
[T8] of Glauber Theory. In the target fragmentation region, on the contrary, the beam 
of independent target nucleons interact with the projectile nucleus, these two results 
coinciding in the central region. The corrections for energy conservation play here a 
very important role if the initial energy is not very high. This approach was used in 
[TO] for the succsessful description of 7r ± , K ± , p, and p produced in Pb-Pb collisions at 
158 GeV per nucleon. 

In the present paper we consider the spectra of secondary baryons and antibaryons, 
as well as their differences, i.e. the net baryon spectra. At low energies, the net baryon 
spectra coincide with the spectra of baryon, while at asymptotically high energies 
they are negligible due to the even signature of the Pomeron trajectory. The energy 
dependence of the net baryon spectra between these two limits strongly depends on 
the production mechanism. 

2.2 String Junction contribution 

In the string models, baryons are considered as configurations consisting of three con- 
nected strings (related to three valence quarks), called String Junction (SJ) [25 ] 126 ] [27 ] 
|2"8"] . this picture leading to some quite general phenomenological predictions. 

The production of a baryon-antibaryon pair in the central region usually occurs via 
SJ-SJ pair production (SJ has upper color indices, whereas anti-SJ (SJ) has lower 
indices), which then combines with sea quarks and sea antiquarks into a BB pair 
[271 [33] . as it is shown in Fig. 2a. 

However, in the processes with incident baryons there exists another possibility to 
produce a secondary baryon in the central region, called SJ diffusion. The quantitative 
description of the baryon number transfer due to SJ diffusion in rapidity space was 
obtained in [TJ and following papers g[ U 1213 I2H 1221 [21] . 

In the QGSM the differences in the spectra of secondary baryons and antibaryons 
appear for processes which present SJ diffusion in rapidity space. These differences 
only vanish rather slowly when the energy increases. 

To obtain the net baryon charge, and according to ref. [19], we consider three 
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Fi gure 2: QGSM diagrams describing secondary baryon production: (a) usual BB central production 
with production of new SJ pair; (b) initial SJ together with two valence quarks and one sea quark; 
(c) initial SJ together with one valence quark and two sea quarks; (d) initial SJ together with three 
sea quarks. 

different possibilities. The first one is the fragmentation of the diquark giving rise to 
a leading baryon (Fig. 2b). A second possibility is to produce a leading meson in the 
first break-up of the string and a baryon in a subsequent break-up [121 El] (Fig. 2c). 
In these two first cases the baryon number transfer is possible only for short distances 
in rapidity. In the third case, shown in Fig. 2d, both initial valence quarks recombine 
with sea antiquarks into mesons, M, while a secondary baryon is formed by the SJ 
together with three sea quarks [5j [191 ES] . 

The fragmentation functions for the secondary baryon B production corresponding 
to the three processes shown in Fig. 2b, 2c, and 2d can be written as follows (see [19] 
for more details): 

G*(z) = a N -vf q -z 2 \ (9) 
G B qs {z) = a N -vf s -z 2 -(l-z), (10) 
Gl(z) = a N -e-vf s -z 1 -^ -{l-z)\ (11) 

where a N is the normalization parameter, and v qq , v qs , vf s are the relative probabilities 
for different baryons production that can be found by simple quark combinatorics 
[351 [37]. These probabilities depend on the strangeness suppression factor S/L, and 
we use S/L = 0.32 following [38] . 
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The contribution of the graph in Fig. 2d has in QGSM a coefficient e which deter- 
mines the small probability for such a baryon number transfer. 

The fraction z of the incident baryon energy carried by the secondary baryon de- 
creases from Fig. 2b to Fig. 2d. Only the processes in Fig. 2d can contribute to the 
inclusive spectra in the central region at high energies if the value of the intercept of the 
SJ exchange Regge-trajectory, asj, is large enough. The analysis in [21] gives a value of 
as j = 0.5 ± 0.1, that is in agreement with the ALICE Collaboration result, asj ~ 0.5 
[55] , obtained at LHC. In the calculations of these effects we use the following values 
of the parameters [24] : 

a s j = 0.5 and e = 0.0757. (12) 



2.3 Contribution from interaction with clusters 

In the case of interaction with a nuclear target some secondaries can be produced in 
the kinematical region forbidden for the interaction with a free nucleon. Such processes 
are called the cumulative ones, the simplest example being the production of secondary 
nucleons in the backward hemisphere in the laboratory frame. 

Usually, the cumulative processes are considered as a result of the coherent inter- 
action of a projectile with a multiquark cluster, i.e. with a group of several nucleons 
which are at short distances from each other that appears as as a fluctuation of the 
nuclear matter [291 HO] ■ 

The inclusive spectra of the secondary hadron h in the central region is determined 
at high energies by double- Pomeron diagrams [4T] . The case of pp collision is shown 
in Fig. 3a. In the case of proton-nucleus collisions two different possibilities exist, the 
interactions with individual target nucleons (Fig. 3b) and the secondary production on 
cluster (Fig. 3c). For nucleus-nucleus collisions there are four possibilities (Figs. 3d, 
3e, 3f, and 3g). The one in Fig. 3g corresponds to a new process where a secondary 
hadron is produced by the interaction of two clusters. 

It was shown in refs. [30, 31] that in the case of secondary production from the 
cluster fragmentation, the inclusive spectra can be calculated in the framework of the 
QGSM with the same quark and diquark fragmentation functions. The only difference 
comes from the quark and diquark distributions, Ug(x,n,k) and Ug l q (x,n, k), where 
k is the number of nucleons in the cluster. The distributions uf(x,n,k) can also be 
calculated by using the Reggeon counting rules, and they take the form: 

u% u (x,n) = C uu ■ x aR ~ 2aB+1 ■ {k - x)~ aRN+mi 
uti(x,n) = C ud -x a »- 2aB . (k - x)- a * +N+m * 
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Figure 3: Reggeon diagrams for the different possibilities corresponding to the inclusive spectra of 
a secondary hadron h produced in the central region in (a) pp, (b,c) pA, and (d— g) AB collisions. 
Pomerons are shown by wavy lines. 
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where 

N = 2 ■ (1 - a B ) ■ (k - 1) . (14) 

The function uf(x,n) does not depend on k. All these functions uf(x,n,k) are nor- 
malized to unity and the values mi and m-i can be found from conditions in Eq. (2). 

This approach was succesfully used in [301 |3T] for the description of cumulative 
particles produced in hA collisions. In the present paper we use it for describing the 
enhancement of strangeness production on nuclear targets in the central region. 

The probability to find a proton in the backward hemisphere in high energy pA 
collisions reach values up to 10%. Keeping in mind that it is only a part of cumulative 
processes, in the numerical calculations and for every hadron-nucleon interaction we 
take the probability to interact with a cluster V d , and the probability to actually 
interact with a nucleon 1 — V cl . We use the numerical values: 

V cl = 0.2 (for mediumA) , V cl = 0.3 (for heavyA) , (15) 

as the maximal reasonable values of V cl . 
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2.4 Inelastic screening (percolation) effects 

The QGSM gives a reasonable description [5J EJ QUI H2] of the inclusive spectra of 
different secondaries produced both in hadron-nucleus and in nucleus-nucleus collisions 
at energies ^snn = 14—30 GeV. 

At RHIC energies the situation drastically changes. The spectra of secondaries 
produced in pp collisions are described rather well [5], but the RHIC experimental data 
for Au+Au collisions [431 HI] give clear evidence of the inclusive density suppression 
effects which reduce by a factor ~0.5 the midrapidity inclusive density, when compared 
to the predictions based on the superposition picture [H2 HS1 HZ]- This reduction 
can be explained by the inelastic screening corrections connected to multipomeron 
interactions [32J (see Fig. 4). 




a b c d 



Figure 4: (a) Multiperipheral ladder corresponding to the inclusive cross section of diagram (b), and 
(c) fusion of several ladders corresponding to the inclusive cross section of diagram (d). 

At energies ^/snn < 30—40 GeV, the inelastic processes are determined by the 
production of one (Fig. 4a) or several (Fig. 4c) multiperipheral ladders, and the corre- 
sponding inclusive cross sections are described by the diagrams of Fig. 4b and Fig. 4d. 

In accordance with the Parton Model [451 145] . the fusion of multiperipheral ladders 
shown in Fig. 4c becomes more and more important with the increase of the energy, 
resulting in the reduction of the inclusive density of secondaries. Such processes cor- 
respond to the enhancement Reggeon diagrams of the type of Fig. 4d, and to even 
more complicate ones. All these diagrams are proportional to the squared longitudinal 
form factors of both colliding nuclei [32J , so their contribution becomes negligible when 
the energy decreases. Following the estimations presented in reference [32], the RHIC 
energies are just of the order of magnitude needed to observe this effect. 

However, all quantitative estimations are model dependent. The numerical weight 
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of the contribution of the multipomeron diagrams is rather unclear due to the many un- 
known vertices in these diagrams. The number of unknown parameters can be reduced 
in some models, and, as an example, in reference [32] the Schwimmer model (50] was 
used for the numerical estimations. Also, in [51] the phenomenological multipomeron 
vertices of eikonal type were introduced for the summation of the enhancement dia- 
gram. 

Another possibility to estimate the contribution of the diagrams with Pomeron in- 
teraction comes [521 E31 EH E51 EE] from Percolation Theory. The percolation approach 
and its previous version, the String Fusion Model [571 ESI EH] , predicted the multiplicity 
suppression seen at RHIC energies, long before any RHIC data were measured. 

New calculations of inclusive densities and multiplicities in percolation theory both 
in pp [60| [6Tj . and in heavy ion collisions [6T[ [62] . are in good agreement with the 
experimental data in a wide energy region. 

In order to account for the percolation (inelastic screening) effects in the QGSM, 
it is technically more simple [23] to consider the maximal number of Pomerons n max 
emitted by one nucleon in the central region that can be cut. These cut Pomerons lead 
to the different final states. Then the contributions of all diagrams with n < n max are 
accounted for as at lower energies. The larger number of Pomerons n > n max can also 
be emitted obeying the unitarity constraint, but due to the fusion in the final state (at 
the quark-gluon string stage), the cut of n > n max Pomerons results in the same final 
state as the cut of n max Pomerons. 

By doing this, all model calculations become very similar to the percolation ap- 
proach. The QGSM fragmentation formalism allows one to calculate the integrated 
over pt spectra of different secondaries as the functions of rapidity and Xp. 

In this frame, we obtain a reasonable agreement with the experimental data on 
the inclusive spectra of secondaries produced in d+Au collisions at RHIC energy [23] 
with a value n max = 13, and in p+Pb collisions at LHC energy [53] with the value 
n max = 23. It has been shown in [M] that the number of strings that can be used for 
the secondary production should increase with the initial energy. 
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3 Numerical results 



3.1 Net baryon spectra 

One example of the QGSM description of the x^-spectra of secondary protons and 
antiprotons measured in pp collisions at 158 GeV/c by NA49 Collaboration [65] is 
presented in Fig. 5. 




Figure 5: QGSM a^-spectra of secondary protons and antiprotons produced in pp collisions at 158 
GeV/c compared to the experimental data by the NA49 Collaboration [65] , 

The QGSM results for net proton (p — p) and net A-hyperon (A — A) productions 
at 200 GeV per nucleon are compared to the experimental data by the NA35 Collab- 
oration [66] in Figs. 6 and 7. 

The data for net baryon production by proton beam interaction with 32 S and Au 
nuclear targets are presented in Fig. 6 in the central and beam fragmentation regions as 
function of rapidity in the laboratory system. The absolute normalization of dn/dy in 
all cases is determined by the data of proton and antiproton production in pp collisions 
at similar energies. 

The results of the QGSM calculations without SJ and cluster contributions are 
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Figure 6: Net proton p — p (upper panels) and net A-hyperon A — A (lower panels) production in 
p— 32 S (left panels) and p— 197 Au (right panels) collisions at 200 GeV per nucleon. Solid curves show 
the QGSM calculations with both SJ and cluster contributions, dashed curves with SJ contributions 
but without the cluster ones, and dotted curves without both SJ and cluster contributions. 

shown in Fig. 6 by dotted lines. Dashed lines show the same calculations with SJ 
contributions but without the cluster ones, and solid lines show the results with both SJ 
and cluster contributions. The corrections for very high energy interactions described 
in Subsection 2.4 are negligible at this energy. 

In the case of net proton production in p— 32 S collisions, all three curves are close 
to each other and they are in reasonable agreement with the experimental data. In the 
case of p— 197 Au collisions the number of net protons is too small at small rapidities, 
what can be explained by the influence of the target fragmentation region. The nuclear 
cluster contribution, which is important mainly in the beam fragmentation region, 
seems to be too large. 

In the case of net A-hyperon production the experimental error bars are rather large, 
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Figure 7: Net proton p—p (upper panels) and net A-hyperon A — A (lower panels) production in 32 S - 
32 S (left panels) and in 32 S - 197 Au (right panels) collisions at 200 GeV per nucleon. Solid curves show 
the QGSM calculations with both SJ and cluster contributions, dashed curves with SJ contributions 
but without the cluster ones, and dotted curves without both SJ and cluster contributions. 

and one can talk of general semiquantitative agreement of the QGSM calculations with 
the data in the lower panels of Fig. 6. 

For the sulphur 32 S beam shown in Fig. 7, the theoretical calculations are in rea- 
sonable agreement with the data for net proton production, whereas the data on net 
A-hyperon production are systematicaly higher than all calculated curves. In any case, 
though, all disagreements are not large than ~30%. 

The data of the NA35 Collaboration [66] are compatible with the results of the 
NA44 Collaboration [67]. 
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3.2 A and A production in midrapidity region 

The NA49 Collaboration obtained experimental data (SB] for yields of A and A hyperons 
in midrapidity region (\y\ < 0.4) in the central C+C, Si+Si, and Pb+Pb collisions (5% 
centrality), at 158 GeV per nucleon. These results are presented in Table 1, together 
with the QGSM results obtained for the same rapidities and centralities. 



(GeV) 


Reaction 
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Experiment dn/dy 


17.2 


C +C ->■ A 


0.237 


0.24 ±0.01 ±0.04, [68] 




C + C -)• A 
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0.064 ± 0.003 ± 0.010 [68J 
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Si + Si ->• A 
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0.68 ±0.04 ±0.13, [68] 




Si + Si ->• A 
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0.16 ±0.007 ±0.038, [SB] 
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Pb + Pb -»■ A 


9.4 


12.9 ±0.7 ±1.5 [68] 




Pb + Pb ->■ A 
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1.4 ±0.3 ±0.2 [68J 


17.2 


p + Be -> A 


0.034 


0.034 ±0.0005 ±0.003 [69J 
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0.011 ±0.0002 ±0.001 [69J 


17.2 


p + Pb -)> A 
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0.060 ± 0.002 ± 0.006 [69J 
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p + Pb A 


0.0019 


0.015 ± 0.001 ± 0.002 [69J 


17.2 


Pb + Pb -> A 


9.4 


18.5 ± 1.1 ± 1.8 [69] 




Pb + Pb ->■ A 


2.05 


2.44 ±0.14 ±0.24 [69J 


62.4 


Au + Au -> A 


11.1 


15.7 ±0.3 ±2.3 [70] 




Au + Au -)> A 


8.2 


8.3 ± 0.2 ± 1.1 [70J 


200 


Cu ± Cu -> A 


3.82 


4.68 ± 0.45 [71] 




Cu + Cu -> A 


3.34 


3.79 ±0.37 [71] 


200 


Au + Au A 


14.2 


14.8 ±2.4 [71] 




Au + Au A 


12.1 


11.7±0.9 P] 


3000 


Pb + Pb -> A 
Pb+Pb -> A 


36.2 
35.6 





Table 1. Experimental NA49 [68], NA57 [69], and STAR jZQl [7T] data for A and 
A production at 158GeV per nucleon, and at STAR energies, and the corresponding 
description by the QGSM. 

On the other hand, the NA57 Collaboration obtained the experimental data [69] for 
A and A yields in midrapidity region \y\ < 0.5 in the minimum bias p±Be and p±Pb 
interactions, and in central (5% centrality) Pb±Pb collisions at 158 GeV per nucleon. 

Unfortunately, the data by the NA49 and NA57 Collaborations are not compatible, 
as one can see from Table 1 , where the values of dn / dy for different hyperons measured 
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by one collaboration are far outside the error bars of the corresponding values published 
by the other collaboration for the same centrality. This is probably due to different 
experimental event selection. . 

Here again one can see that the calculated yields of A and A are in agreement with 
experimental data on the level of 20—30% accuracy. 

Hyperon production at higher energies in midrapidity region was also measured 
at RHIC. The data by the STAR Collaboration pU [71] for Au + Au and Cu+Cu 
collisions at ^s^n = 62.4 GeV and 200 Gev are presented in Table 1. We also give in 
Table 1 the QGSM predictions for central Pb+Pb collisions at the LHC energy ^snn 
= 3 TeV. 

4 Conclusion 

The QGSM provides a reasonable description of nucleon and A, as well as their an- 
tiparticles, production in nucleon-nucleus and nucleus-nucleus collisions at high ener- 
gies. The level of numerical accuracy is of about 20—30%. Part of the uncertainty is 
connected to discrepancies among the different experimental data. 

Inclusive densities of H+ hyperons are reasonably reproduced for the cases of p+Be 
and p+Pb collisions (69], but they are several times underestimated in the case of 
central Pb+Pb interactions [69J. For fi+ production in central Pb+Pb collisions the 
disagreement is larger than one order of magnitude. The physical reasons for these 
observed disagreements will be discussed in a separate paper. 
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